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ABSTRACT 


The transfer of kinetic energy by individual wave 
numbers is computed at various latitudes using data ob- 
tained from a harmonic analysis of the hemispheric 500 
millibar pressure height. Linear correlations are sought 
between the computed values of kinetic energy transfer and 
computed 24-hour changes of kinetic energy of the pertur- 
bations for the same and for adjacent wave numbers. 

The writer wishes to express sincerest appreciation 
to Professor Frank L. Martin for his suggestion of the 


topic and for his guidance during the study. 
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1. DGetxneducti on. 

One phase in the study of the general circulation of 
the atmosphere has been the investigation of mechanisms by 
which the kinetic energy balance is maintained in spite of 
continual frictional dissipation Widger [14], Kao [6], 
and others have established that a mean northward trans- 
port of kinetic energy in mid-latitudes is effected by the 
lange scabe: perturbations in the atmesphere. Saltzman 
[9, 10] has investigated the spectral distribution of per- 
turbation kinetic energy. Saltzman and Fleisher [1l] have 
Studied the rate of transfer of kinetic energy between the 
perturbations and the mean flow on a daily basis, while 
Benton and Kahn [1] have conducted a Similar investigation 
on a Seasonal basis. Saltzman and Fleisher [12] have also 
estimated the interchange of eddy kinetic energy between 
groups of perturbations, assuming that the total kinetic 
energy interchange is zero. 

It is the purpose of this study to correlate the hori- 
Zental geostrophic transfer of kinetic energy by wave num- 
ber with 24-hour changes of kinetic energy of the same and 
different wave numbers existing in the same and nearby 
latitude belts. The data used are from a Fourier analysis 
of 500-mb contours at successive 5-degree latitude belts. 
This Fourier analysis, which covers the winter season of 
1947-48, was made available by the U. S. Navy Weather 
ResearcheFacpilaty;, Norfolk. ~Further discussion concerning 


the data appears ina later section. 





2. The energy transfer function. 

Using the equations of quasi-horizontal motion, 
Saltzman [9] has developed an expression for the change of 
eddy kinetic energy within a closed mass of air, /‘\, eggs 


the entire atmosphere, which is 


~ / = 7 Googe | to ee 
2 | KE aoe uv Q 2 (Sa ee ET 
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Here D is latitude, A longitude, pressure, and 2, the 


earth-radius. The eastward and northward components of the 
wind are WU and V respectively, while the vertical motion 
is W. Primes are associated with perturbation quantities 
and bars with averages around a latitude circle. Thus, for 


example, 


= 
= &[0)% 
3 

indicates the barred-type average. The mass element dm in 
equation (2.1) will be discussed in a later section. 
Saltzman also shows that the change in the mean kinetic 
energy, & (KE Sm , contains the first integral on the 
right in (2.1), but with the opposite sign affixed. 

It may be seen that the first integral in equation 
(2.1) contains terms which are products of the macroscopic 


Reynolds stresses and the derivatives of the mean flow. 





The remaining integrals represent contributions to the eddy 
kinetic energy by conversion of "eddy available potential 
energy" and frictional dissipation, respectively. Inasmuch 
as this study is concerned only with the effects of hori- 
zontal transfer of kinetic energy, no further discussion of 
the geopotential or frictional terms will be entertained. 
If the geostrophic approximation is made, as is done 
in this study, the mean meridional component “Viewi We Vemmiesh. 
Furthermore, since it is impossible to determine vertical 
motions from the Fourier analysis of the height field at a 
Single level, those terms involving GC)” will be neglected. 
With these simplifications, the eddy kinetic energy 


change may be expressed in the form 


2 i ier VV Ie CIOS 2 (u.) 
7 K Sa Lene Uae oo ém (2.2) 
e M 


For an interpretation of the integral appearing in 
equation (2.2), a vertical element of area, dS = a cosh 6102, 
contained in a latitude wall, will first be considered. 


Here the eastward momentum per unit mass of air is 


Any transport of eastward momentum normal to 68 must be 


brought about by the meridional wind component 
aos nee 
Vs= VT V 


With the geostrophic approximation, the longitudinally- 
averaged meridional transport of eastward momentum due to 
large scale turbulence is given, after Haltiner and 
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Mamtin [4], 
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The integral in equation (2.2) will be denoted ia 
and will be termed the "eddy transfer function." It will 
eventually involve contributions to the kinetic energy 
change by different wave numbers. A positive value for 
ae represents an exchange of kinetic energy from the mean 
flow into that of the eddy flow. 

3. The energy-transfer spectrum. 

Since U and V may be considered real and continuous 
functions of longitude ~ , each having unique values at all 
points of a given latitude circle, they may be expanded in- 


to Fourier series of the form 


Ula) = a0) + > Qyh) cos nr + Din) Sin na (3.1) 
where = 
a,(o) = cf uo 52 
oO 
2 | . 
_ ‘sc hHaw See 
q, (h) = Tr | A A 
° on 
b(n) = #[ uc 
Likewise 
Via) = GO) + DA, In) cos na + Din) sin na 3 3) 


where a_(o) p WARE oy, and b(n) are defined similarly to 


equations (3.2). ™quations (3.1) and (3.3) may also be 
4 








Gapressed, in tésms Of amplitude and phase angle, as 
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with the amplitude and phase angle of equation (3.4a) given 


by 
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Now the copower series of U and V may be defined, 


after Kahn [5], 


oh if uaVadsa = a0) ao +> Com (5.6) 


where it may be verified upon combining equations (3.4) 
and (3.6), that the copower spectrum, or cospectrum, has 


the form 


Coin = eHmHim cos nam) ~ &)) (3.7) 


Here the symbols H(n) , and oye are defined as in equa- 
tions (3.5). It may also be noted that for the latitudinal- 
mean values of Uta) or vr (a) , the cospectrum Ge in 

equation (3.6), becomes simply the power spectrum of ua) 


A 


or Vv (a) 


fs £ 
Rn) = +H (n) or a thi = 41H in) (3.8) 
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Now the value of OF ae any given wave number may be shown 


to be independent of its values at any other wave number. 





Equations (3.7) and (3.8) may therefore be used to ob- 
tain the contribution of a specific wave number to 
either the cospectrum or the power spectrum, as the case 
may be. 

A useful expression for ava may be obtained by 


subtracting Q (0) {0 from the left side of equation 


(3.6), after which 


uv’ = sae (h) 


Making this substitution into the transfer equation (2.2) 


gives 
i. - (5c (ny $282 * (caso) a $ Tm (3.9) 
M he | nat 


Thus a may be expressed as a sum of contributions for 
each wave number. 

A more convenient form of equation (3.9) will now be 
derived. Making use of the geostrophic approximation to 
the horizontal flow, the eastward and northward compo- 


nents are 


CF 


J 
fa oo 
(3.10) 
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Due 
facosSd oa 


Here 9 is the acceleration of gravity, if the coriolis 
parameter, and Z contour height. The other symbols are 


as previously defined. Denoting the n-th harmonic of 


6 





Bie Momeme fre 1d “at lLatutude D as 


Zin) = Ath) Cos n(n - Yn) Rod 


where An) and (n) are the amplitude and phase, respec- 


tively, at latitude dD , equations (3.10) are expressible 


as 
, 
din) = = Awa cos nla- te) + NARA rn lao) 3.12) 
| w 
Vin) = elite eee nla- Ya) 3.13, 
~ $acos 6 


Henceforth the wave number !) will be understood in Aw) ana 


Yn). Equation (3.12) may also be written in the form 


a= -2) 4/64) 


where the angle Voadus given by 


ae (25 / nh ) 3.15) 


and is therefore the phase difference between Ug, and 


Vo . Using the expressions (3.13), (3.14), (3.15) and 
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(3.7), equation (3.9) takes the form 
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Boaucieion (I2l6)ewill be termed the energy transfer spece- 
trum. By simplifying the quantity within the square 
brackets, the contribution for a single wave number 


becomes 
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4. The 24-hour change of eddy kinetic energy. 

An approximation to the change of eddy kinetic energy 
in a closed system resulting from all of the processes of 
equation (2.1) is readily made. The eddy kinetic energy 


per unit mass at any position (ao) is given by 


g 2 2 
K E(n) = ¥en) —- (uin_4 'in) (4,1) 


Using equations (3.8), (3.12) and (3.13), one obtains the 
latitudinally-averaged or spectral values at latitude dD 
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Substitution of equations (4.2) into (4.1) yields the 
zonally-averaged geostrophic eddy kinetic energy per unit 


mass. For a closed system of mass jae, 
- g {1 Trane soe, inl FI _ 
KCl = ee alco 14 i (235 6) | m 4.3) 
M 
The mean 24-hour change is found uSing the simple differ- 


ence, shown symbolically, 


A KEW oe KEM aaete KE (hy (4.4) 
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5. Finite-difference forms. 


The element of mass 5m of equations (3.17) and (4.3) 


may be expressed initially as 


U1 
~~ 


dm = parcosh 25d? (5. 





2s = ae 
— & - - - eee a el 
ti ess a 


which, upon substitution from the hydrostatic equation, 


becomes 





se — $a So Sp (5.2) 


where bp, the pressure increment is considered positive. 

In a series of papers, Saltzman [10], and Saltzman and 
Fleisher [1ll, 12] have been careful to simulate a closed 
system by uSing a large portion of the northern hemisphere 
at 500 mb, for example, from 15° N to 80° N. In this thesis, 
it seemed advantageous, from the point of view of a possible 
forecast tool, to investigate kinetic-energy changes in ten- 
degree latitude belts, even though such belts may not con- 
stitute closed systems. 

The major purpose of this study is to investigate pos- 
Sible relationships between a we and OKEM at the 500-mb 


level. Combining equations (3.17) and (5.2) leads to 
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ou ra AN ose Sa) faze) 
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Integration of (5.3) with respect to longitude will increase 
the constant term by a factor of 2/] . The remaining inte- 


grations are with respect to pressure, where 


centered at 500 mb, and with respect to latitude, where 
2A = 10 degrees, or 7| /18 radians. 
Employing finite differences and the geostrophic 


approximation, the quantity & in equation (5.3) becomes 
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is the height of the mean height contour at the 


latitude indicated by the subscripts. Moreover LU, at the 


central latitude D , may be expressed as the arithmetic 
mean 
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Substitubing (5.4) and (5.5) into the integrated foxzm of 


equation (5.3), the. finite difference form may be stated as 
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Combining equations (4.3) and (5.2) yields 
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which, when integrated as before, takes the finite-difference 


form 
KE a Pos Ae - 3) n 2A (tases Pee) 
‘EWU ~ “ag | sino 2A d —— 
hide Vi 20 v (58) 
cosh |} 
where 


A=4(A,,. 7 cA\y + A o ee aad AY te 


It should be noted that the same integrations have been 


performed in arriving at equation (5.8) as for (5.6). The 
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UnLES Of K E(n) are™théreforne those of Kinetic Gnergy per 
mb and per 10° latitude belt. 
6. Data processing and computation. 

Data for the evaluation of equations (5.6) and (5.8) 
were taken from a Fourier analysis of the 500-mb pressure 
Surface which was made by the Naval Weather Research 
Facility [13]. The Facility used the Historical Series 
500-mb hemispheric charts for the winter season of 1947- 
48. Mean 500-mb heights, and the amplitude and phase for 
each of the first 18 wave numbers were computed Values 
were available for five-degree intervals from 15 degrees 
to 70 degrees North latitude. Data selected for this 
study were for the 33-day period from 31 December through 
1 February. 

Because of time limitations, only waves 1 through 12 
were considered. However, results of Duggan [3] and 
Melhorn and Le Dew [7] indicate that contributions from 
wave numbers higher than ten may be expected to be rela- 
tively small. 

In the case of phase-angle data, further refinement 
was necessary. Phase angles were given in terms of dis- 
tance from the prime meridian to the ridge line of the 
wave in question. Because of the method of tabulation of 
phase angles used in the original harmonic analysis, this 
value did not always pertain to the ridge nearest the 
meridian. Thus, if the phase difference obtained from the 


reported values was in excess of one-half wave length, it 
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was necessary to subtract one wave length from the orig- 
inal difference. Fig. 1 schematically depicts this process 
Gf identification. 

Computation of equations (5.6) and (5.8) was performed 
for each of five ten-degree latitude belts, centered re- 
spectively at 65, 55, 45, 35, and 25 degrees. The Control 
Data Corporation CDC 1604 computer was employed in these 
evaluations, being ideally suited because of its high speed 
and large capacity. 

7. Results and conclusions. 

It was indicated previously that if integration of the 
transfer equation (2.2) is performed over the entire atmos- 
paene, ANnJeapproxinmete measure of @he transfer of kinetic 
energy between the eddies and the mean zonal wind is ob- 
tagned. JjFrom the evaluation of sthis integral over a 10- 
degree latitude belt, as is done in this study, only a 
palelaleweterEtCation Of~this approximation can be inferred. 
Fig. 2 shows the 33-day mean-energy transfer spectrum 
Tin) at the central latitude of each 10 degree latitude 
belt. A positive value of ai einciasees a contribution £o 
the kinetic energy of the eddy flow from the mean wind, 
whereas a negative value indicates that the eddies supply 
kinetic energy to the mean wind. 

Perhaps the most striking feature of the spectra 
shown in Fig. 2 is the change in sign of the transfer 
between 35°N and 45°N. At the lower latitudes, transfer 
from the eddies to the mean flow is indicated. This is 


in keeping with the findings of Kao [6] who shows a maximum 
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convergence of northward transport of kinetic energy at 
500 mb up to approximately 36°N, which is in the vicinity 
of the zone of maximum westerlies. Thus, at latitude 

25°N and 35 °N the transfer spectrum appears to assist in 
the maintenance of mean zonal winds. It should be noted 
that waves 2 and 3 are the largest individual contributors 
at these latitudes. 

At A5°N, a change in the sense of energy transfer may 
be seen. All waves, with the exceptions of 7 and 11, now 
receive energy from the mean flow. The transfer from mean 
to eddy flow occurs at 55°N and 65°N, although at 65°N the 
magnitude of the exchange is greatly reduced. 

dimes spectral estimate for the full latitude range 20°N 
to 70°N is shown in Fig. 3. In this instance the sign of 
the energy transfer spectrum has been changed to - | (| 
to facilitate comparison with the results of Saltzman [10]. 
In Fig. 3, a positive sign now indicates a contribution to 
the mean wind by the eddies. 

Results of this study agree roughly with those of 
Saltzman [10]. The major difference between the two 
spectra lies in the net transfer for all of the waves. 


Saltzman's results show an overall transfer from the eddies 
si el a F | > CL, whisle ain the 


to the mean flow, | Lz 
: Niney / 


present study there is a net transfer of kinetic energy 


+ 


T 


from the mean to the eddy flow, |, / | <, AL. This 
is not unexpected in the light of Saltzman and Fleisher's 
Peaimauoated Values of the standard deviation of the transfer, 


which, in general, exceed the value of the spectral estimate. 
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It is reasonable to suppose that changes in Tin at 
a given wave number may be accompanied by changes in the 
kinetic energy associated with other wave numbers at the 
same latitude and at adjacent latitudes. Accordingly, 
values of Tiny at latitude 45 degrees were linearly cor- 
related with those of 7, .. | for wave numbers n = l, 


c | XK 
2s « sas In addition, Ga) was correlated with values of 


LN} at adjacent wave numbers nt 1 andn* 2. Corre- 
sponding values of /\ at latitudes 35 and 55 degrees 
were also correlated with the value of i at 45 degrees. 
faowe eves marsehematic representation of ths correa-— 


tion system. 


hat 55° WKE(n-2) asKE(n-1) «4KE(n) 4aKE(n 1) KE(n 2) 

meeag? =aeHi(in- 2) KE(n-1) 2KE(n) -KE(n 1) «<KE(n 2) 

Lat 35° AKE(n-2) AKE(n-1) AKE(n) sKE(n 1) AKE(n 2) 
Table 1. 15-way correlation scheme. 


T*(n) is correlated with each 
AK.E. shown. 


Table 2 shows the results of the 15-way correlation 


rc > * 


system in which | (nv at 45 degrees is correlated with each 
of the 15 sets of values of AKE(n,n‘l,n=2) at latitude 35, 


45, and 55 degrees. The values of AKE are centered 12 


Weurs™prewor to the time of "“'™™. Tables 3, 4, and 5 Show 
the results with \'.[ centered 12, 36, and 60 hours respec- 
tively after the time of | inh, The tables are arranged in 


the schematic fashion shown in Table l. 
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eupexcellent statistical routéimey written for the 
CDC 1604 computer by Campbell [2], was employed in the 
computation of the correlations. 

In proceeding to the correlation analysis, it was 
meGessary £O Establish confidence limits for the correla- 
tion coefficient. The number of independent pairs in all 
the correlations was 30. According to the z-test, as : 
given by Panofsky and Brier [8], a correlation coefficient 
of .30 is Significant at the 90% level of belief. 

The physical principle to be tested is the conversion 
of the transfer-energy by wave number n into the kinetic 
energy of wave number n, nil, ni2. Saltzman [9, equa- 
tion (49)] has derived a transform of equation (2. ljwe His 
result includes all of the terms which appear in (2.1), in 
spectral form, and in addition, terms representing inter- 
action of kinetic energy from wave number n to all other 
waves. Ina later paper, Saltzman and Fleisher [12] have 
estimated the rate of energy interaction at 500 mb. They 
have shown that energy transitions tend to occur between 
a middle group of waves, (n = 6 through 10) called "cy- 
clone waves," into the longer wave group (n = 1 through 5) 
and the short wave group (n = 11 through 15). The present 
study is more limited in scope, in that the range of waves 
considered is only (n - 2) through (n + 2). On the other 
hand, investigation by means of correlations between speci- 
fic wave numbers may make available more detailed criteria 
for the energy exchange between waves at latitudes Bick? 


“A 


5 and 56/0 in terms of at 45 degrees latitude. 
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For the purpose of such an investigation correlation 
coefficients from a particular row and column of each of 
Tables 2, 3, 4, and 5 were selected in accordance with the 
criteria that each set of four successive coefficients 
contain at least: 

(a) two significant correlations 
i (b) one significant correlation together with 
a 24-hour or 48-hour difference in corre- 
tatecneor = .60. 
All together there were 47 sets of successive correlations 
in Tables 2 through 5 which satisfied criteria (a) or (b). 

As an example, if the row and column just mentioned 
correspond to (ll, 9), each of the four tables would be 
entered with the central n-value, n = ll, and column 
n- 2:29. This particular correlation coefficient would 
be symbolized ¥ .. . If a row corresponding to the higher 
latitude (55°N) is selected, a bar is placed over the | 
(e.g., Y ), whereas if the lowest latitude is used, a bar 
is placed below the Y (e.g., ~ ). For a choice of a 
central latitude, no bar is employed. 

For example, J \, 9 signifies the correlation between 
nee: computed at 45°N versus AKE(9) computed at latitude 
35 degrees. Of this particular set of four correlations, 
three were significant at the 90% level or higher. The set 
of four is plotted in Fig. 4a. are shows a similar 
pattern, but in this case has a significant positive cor- 


relation at time t -— 12 hours. This initially positive 


Significant correlation tends to confirm the hypothesis 
——y) 


that an above-normal value of /! (N) is associated with an 
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above-normal value in MAKL(n) at the same wave number. 
At time t + 12 hours, however, the value is negative, os- 
cillating to positive and negative again on the succeed- 
ing two days. These results are consistent for all cor- 
Hediatieonmses meetbingsthescmitesija ofa) o 4b) an mp O 
and suggest that above-normal value of Tin) tends to be 
followed by below-normal change in perturbation kinetic 
energy twelve hours later. Several of the correlation 
sets examined in this group displayed this two-day peri- 
odicity. Certain of the other sets in the short-wave 
group lose significance at time t + 36 and t + 60 hours 
even though they may have met the criteria (a) or (b). 

At the opposite end of the wave number scale, there 
is again the tendency for AKE (2 or 3) to be above-normal 
with a high value ot | fel. as shown in Fig. 4b. Again 
AKE (2 or 3) becomes negative, but in this instance, two 
days later, indicating the periodicity of the kinetic 
energy changes at waves 2 and 3 is apparently of the order 
of 4 days. 

At wave number 5, there is a marked tendency for a 

WE. (5) to be negative at t + 12 hours. A tendency for 

Z “rT (6) to be negative at this same time also exists as 
Shown in Fig. 4c. This suggests a positive correlation 
between the two. However, the exact nature of any such 
Pebatronsiip LS difficult to determine by the technique 
emproyed here. 


At intermediate wave numbers, as represented by n = 8, 


tends to be positive at t + 12 hours. At the same 


Zk 


> 


C2 CC 
one 





time, however, 3 6 , and 1 are of opposite sign as 
may be seen in Fig. 4d. This result points to a general 
agreement with the results of Saltzman and Fleisher [12], 
namely, when a decrease of kinetic energy occurs in the 
cyclone waves, kinetic energy tends to be exchanged both 
upward and downward in the wave number scale. A schematic 


illustration of this exchanged energy, after Saltzman and 


Plersm@er [l2)Ltseshownwine Fig. 5. 





Figure 5. Estimate of the kinetic jenecmgy 
exchange between cyclone waves and 
longer and shorter waves. After 
Saltzman and Fleisher [12]. 


A direct correlation analysis of the kinetic energy 
changes between the various waves would have yielded seven 
times aS many correlations as the present scheme, so that 
analysis of the results would be extremely involved. Also 
the results would have had no relationship to T* (n), which 
function affords the possibility of being a predictor 
based on physical considerations. 

It is necessary to comment on the relatively small 


number of Significant correlation coefficients. For 
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Gxempre, in’ Tables 2 and 3, “there occurred 27 and 30 sig- 
nificant values, respectively, where 16 (or 10%) could 
have arisen due to chance. Some of the reasons that a 
more conclusive number of Sigmaficant correlations did not 
appear may be the following: 
(a) Of the several macroscopic eddy stress 
terms in equation (2.1) only the T* term 
was considered. 
(b) Twenty-four hour changes of KE may not 
behave linearly when correlated with Ti. 
(c) In many cases the initial correlation be- 
tween T*(n) and AKE(n) at t - 12 hours was 
negative but became Significantly positive a 
day later. This 1S admdveaqeiye oLerene in 
portance of the process of exchange of kin- 
etic energy between the waves in the manner 
described by Saltzman [12]. Such exchanges 
undoubtedly contribute to the initial state 
of the kinetic energy spectrum. 
(d) There is some “noise" present in all 
Fourier analyses of meteorological parameters. 
There may also have been some undetectable 
reproduction errors in the copy of the 
original data used in this Study. 
In spite of the shortcomings of the linear correla- 
elon technique employed in this study, 2t 1s felesinae ase- 


laemng- sigmaticant” correlation sequences of the type 
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Shown in Fig. 4, led to several valuable results. Among 
these~arre : 
(a) the short period of kinetic energy 
changes at the high wave numbers; 
(b) the longer period changes at lower 
wave numbers; 
(c) the tendency of wave number 8 (in the 
center of the cyclone wave group) to ex- 
change kinetic energy with both waves 6 


and 10. 
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APPENDIX 


The following tables list the results of correlating 
the energy transfer function T* (n) computed at A5°N, with 
the 24-hour Kinetic energy change AKE computed at 35°, 
aS, and 55°N. In addition, AKE is varied with respect to 
wave number and time as indicated by the respective table 
Gaperuems.) Signiticant Correlation coeffirCiénts age de- 


noted by an asterisk. 
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‘ Correlations of ui (n) rat '45ON. and time Ct, Wiea 
mm coe ng” nyt), ears pee «+2, at time t.+ ‘12 hours. 
BRA EA ane. values are “marked ween. an asterisk. 
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Table 4: Correlations of T”(n) at 45°N ard time t with 
:. SRE ae nwt iil, and n 2,ac tame t % 36 hours. 
Significant values are marked with an asterisk. 
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Table 5. Correlations of (7) at 45°Nn and time t, with 
4KE'at n, ni-1l, and n+ 2,at tire t + 60 hours 
Significant values are marked with an asterisk. 

31 























